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ABSTRACT 

Both  analytical  and  experimental  studies  were  conducted  on  a  pro- 
pagating crack  in  a  viscoelastic  material.  By  applying  an  extension  of 
the  "correspondence  principle"  the  stress  and  displacement  at  the  crack 
tip  were  found  as  functions  of  the  crack  tip  stress  intensity  factor. 
An  energy  balance  was  made  for  the  specimen  which  allowed  this  factor 
to  be  related  to  geometric  and  load  parameters.  Fracture  characteriza- 
tion was  then  performed  by  experimentally  relating  the  crack  tip  stress 
intensity  factor  to  the  crack  velocity. 

The  theory  was  applied  to  solid  propellant  fracture  tests;  however, 
the  dewetting  behavior  of  the  propellant  was  not  accounted  for  by  the 
theory. 
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I.   INTRODUCTION 

In   the  elastic  crack  propagation  problem,  it  is  useful  to  define  a 
crack  tip  stress  intensity  factor,  K.   This  factor,  which  includes  the 
effects  of  specimen  loading  and  geometry,  gives  a  measure  of  the  stress 
concentration  at  the  crack  tip.  As  the  load  on  a  specimen  containing  a 
stationary  crack  is  increased,  K  increases  until  a  critical  value  is 
reached,  whereupon  the  crack  rapidly  propagates.  This  critical  value, 
K  ,  is  experimentally  measured  and  found  to  be  a  material  property. 
A  material  may  then  be  characterized  by  this  constant,  for  if  any  com- 
bination of  geometry  and  loading  should  produce  K  ,  crack  propagation 
will  occur. 

For  cracks  in  viscoelastic  materials,  Knauss  and  Mueller  [Ref.  2] 

stiowed  that  in  effect  there  is  a  certain  minimum  loading  which  will 

initiate  crack  propagation;  that  is,  there  is  a  loading  below  which  the 

crack  velocity  is  small  enough  to  be  considered  equal  to  zero.  However, 

the  velocity  of  propagation  will  vary  over  a  range  of  several  orders  of 

magnitude,  increasing  with  increasing  load.  Thus  the  question  to  be 

answered  in  the  viscoelastic  fracture  problem  is  not  simply  will  the 

crack  propagate,  but  rather  at  what  velocity  will  it  propagate.  Knauss 

[Ref.  3]  proposed  that  for  a  given  viscoelastic  material,  K(c)  is  a 

unique  function,  where  c  refers  to  the  crack  velocity.  Therefore  it  is 

natural  to  experimentally  characterize  viscoelastic  materials  by  this 

relation. 

The  material  property,  Kc,.is  unique  for  most  materials  and  geo- 
metries. However,  there  is  a  variation  of  Kc  with  initial  crack  length 
(usually  small  enough  to  be  ignored)  explained  by  consideration  of  the 
growth  resistance  curve  as  shown  in  Ref.  1. 


A  common  method  of  obtaining  viscoelastic  solutions  is  to  use  the 
"correspondence  principle."  Lee  [Ref.  4]  showed  that  viscoelastic 
boundary  value  problems  possess  the  property  that  an  appropriate  elastic 
solution  in  the  Laplace  transform  domain  may  be  inverted  into  the  time 
domain  to  give  the  viscoelastic  solution.  However,  this  "correspondence 
principle"  does  not  in  general  apply  to  problems  where  the  boundary 
conditions  are  mixed  and  apply  over  regions  which  change  with  time. 

Now  looking  at  the  problem  of  the  propagating  viscoelastic  crack, 
note  that  the  crack  edge  boundary  condition  immediately  behind  the  crack 
tip  is  zero  stress  while  the  boundary  condition  on  the  extended  line  of 
propagation  of  the  crack,  immediately  in  front  of  the  tip,  is  zero 
displacement.  Propagation  of  the  crack  implies  that  portions  of  the 
boundary  under  zero  displacement  are  continually  changing  into  zero 
stress  boundaries.  Thus,  the  restriction  on  the  "correspondence 
principle"  mentioned,  above  prohibits  its  application  to  this  problem. 

A  different  solution  method  has  been  used  by  Dietmann,  Knauss  and 
Mueller  [Refs.  2,5,6,7]  in  their  various  studies  of  propagating  cracks 
in  viscoelastic  materials.   In  all  of  these  cases,  the  actual  problem 
was  approximated  by  one  in  which  crack  propagation  was  a  step-wise 
process.  The  resulting  solutions  in  all  cases  contained  the  length 
of  the  crack  increment  used  in  the  step-wise  process  as  a  material 
property;  however,  this  length  does  not  physically  exist  in  the  actual 
problem  and  therefore  seems  to  be  only  the  result  of  the  solution 
technique. 

In  the  solution  of  the  viscoelastic  fracture  problem  developed  in 
this  paper,  the  "extended  correspondence  principle"  will  be  used.  This 
extension  is  due  to  Graham  [Ref.  8]  who  extended  the  classic 


"correspondence  principle"  to  cover  several  cases  where  mixed  boundary 
conditions  applied  over  regions  which  changed  with  time.  Using  this 
extension,  both  the  viscoelastic  stress  and  displacement  in  the 
vicinity  of  the  tip  of  a  continuously  propagating  crack  will  be 
obtained  as  functions  of  the  stress  intensity  factor,  K.  The  remainder 
of  the  solution  will  consist  of  applying  an  energy  balance  to  the 
specimen  and  relating  K  to  the  loading  and  geometric  parameters  of  the 
specimen. 


II.  THEORY 

Consider  a  thin,  narrow  sheet  of  isotropic,  homogeneous,  linear 
viscoelastic  material  loaded  uniaxially  by  aa(X,t)  (where  X  refers  to 
geometric  coordinates  and  t  to  time)  containing  a  continuously  pro- 
pagating center  crack  of  length  2c(t)  as  shown  in  Figure  1. 
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Figure  1.  Specimen  Geometry  and  Loading. 

Assuming  that  the  material  is  thermorheologically  simple,  the  analy- 
sis for  the  stress  intensity  factor  is  made  at  a  constant  reference 
temperature,  T  .  This  solution  will  then  be  valid  at  any  other 
temperature  when  appropriately  shifted. 


A.  CRACK  TIP  STRESS  AND  DISPLACEMENT 

In  this  analysis  a  knowledge  of  the  work  produced  at  the  crack  tip 
will  be  necessary.  This  further  requires  an  evaluation  of  the  stress 
and  displacement  fields  in  the  vicinity  of  the  moving  crack  tip,  which 
is  now  undertaken. 


If  the  sheet  and  loading  are  doubly  symmetric,  a  single  quadrant 
may  be  considered,  as  shown  in  Figure  2,  where  the  lower  right  quadrant 
is  removed  and  replaced  by  the  stresses  which  act  along  the  sheet's 
horizontal  centerline. 


*•  x 


Figure  2.  Single  Quadrant  of  the  Doubly  Symmetric  Specimen. 
Let  the  entire  boundary  of  the  sheet  be  defined  as  8  with  BoU) 
representing  the  uncracked  portion  of  the  x-axis  and  B-i(t)  representing 
the  remaining  boundary.  Then  8  =  6-.(t)  u  S2(t).  The  boundary  conditions 

are: 

0 

0      on  crack  edge 
on  x  =  0,1 


°s(x,t) 


on  8 


on(X,t) 


un(X,t) 


a  (X,t)  on  y  =  b 

a 

0 


on  8-j(t) 


on  82(t) 


(1) 


where  u  refers  to  displacement  and  the  subscripts  s  and  n  respectively 
refer  to  norr.al  and  shear  components.  The  boundary  conditions  of 


Equation  (1)  meet  the  criteria  of  Graham's  "extended  correspondence 
principle." 

However,  to  apply  this  correspondence,  the  elastic  stresses  and 
displacements  in  the  vicinity  of  a  stationary  crack  tip  are  needed. 
They  have  been  found  by  Irwin  [Ref.  9]  to  be: 


o  =  cos  s- 

n  r~o — ^    *■ 


,     .8-36 
1  +  sin  j  sin  y~ 


{ 


r     .  e 
o —  sin  o- 


2  -  2v  -  cos2  | 


(2) 
(3) 


where  the  subscript  1  on  the  stress  intensity  factor  indicates  that  the 
crack  is  propagating  in  the  opening  mode  as  defined  by  Ref.  10. 

Setting  6  =  0  in  Equation  (2)  and  6  =  u  in  Equation  (3)  gives  the 
stress  and  displacement,  close  to  the  tip,  on  B2(t)  and  $-,(t) 
respectively.   It  is  reasonable  to  assume  that  Equations  (2)  and  (3) 
apply  to  slowly  propagating  cracks  where  the  inertia  effects  are  neglig- 
ibly small  in  the  equilibrium  equation.  Changing  to  an  inertial 
coordinate  system,  where  x  =  c(t)  +  x'  (refer  to  Figure  2)  and  assuming 

the  incompessible  case  (v  =  j) ,   the  elastic  stress  and  displacement  in 

2 

the  vicinity  of  a  slowly  moving  crack  tip  become  the  following: 

K,(t) 


oZ   (x,t)  =  - 


h  (t  -  t) 


on  p2(t)  (4) 


on 


u^  (X,t)  =(2  V^  \    K-|(t)  <Jc(t)-x  h(t-t)  crack 

edge 


u^  (X,t) 
Uo  (X,t) 


on  x  =  0,1  ^  e-,(t)    (5) 
on  y  =  b 


o 

In  elastic  theory,  the  stress  intensity  factor  in  these  equations 
is  the  critical  one,  iC|c,  since  the  crack  is  moving.  However,  in  the 
viscoelastic  case,  with  continuous  propagation,  this  definition  of  the 
critical  value  loses  its  significance. 
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where  superscript  e  refers  to  elastic  solutions,  u-,(X,t)  and  u?(X,t)  are 
elastic  displacements  which  are  functions  of  geometry  and  modulus,  h(t) 
is  the  unit  step  function  and  t  is  the  time  at  which  the  crack  tip 
passed  a  point  on  the  horizontal  axis  located  at  x. 

Graham  has  shown  that  Equations  (4)  and  (5)  are  the  solutions  to  the 
viscoelastic  problem  in  the  Laplace  transform  plane  if  the  elastic 
modulus  is  replaced  by  the  transform  modulus,  E->pE  ,  (p).  Therefore, 
the  field  quantities  of  interest  are: 


on    (X,p)  = 


h  (t  -  t) 


on  B2(t)    (6) 


Un  <X'P>  =  2VT^T^ 


where:         /Vptf(t)dt  =   [f(t)] 


/2tt  (x-c(t)) 

^(t)  Jc(t)  -  x  h(t  -  t) 

on  6-,(t) 


(7) 


Since  pD cr(p)  =  VpEr  -j(p),  where  D   refers  to  creep  compliance, 
the  viscoelastic  solution  inverts  back  into  the  time  domain  with  the 
stress  distribution  identical  to  the  elastic  case  and  the  displacement 
at  the  crack  tip  edge  given  by: 


un  (x't]  =  2V~F   Dcr(t)  Ki(°Wc(0)  '  x  h(0"^ 


"cr*1  "T»  57 


^U)\cU)   -  x  h(x  -  t) 


J 


,(8) 


Since  h(0  -  t)  =  0  and  h(x  -  t)  =  0  for  x  <  t,  Equation  (8)  may  be 
rewritten  as: 


un   b-V-zfirjXrlf^k    ¥T>V^ 


dx 


(9) 
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Time  t  is  seen  to  be  equal  to  time  t  minus  the  time  which  it  took 
for  the  crack  to  run  between  x  and  c(t),  that  is: 


■:  *=  t  -  c^  "  x  (10) 


m 

where  c  is  the  mean  velocity  in  the  interval  [x,c(t)].  Thus,  time  t 
accounts  for  the  history  which  has  affected  the  displacement  at  point  x. 
Or  to  put  it  another  way,  the  crack  displacement  at  point  x,  at  the 
present  time  t,  has  been  creeping  open  ever  since  the  crack  tip  passed 
that  point  at  time  t. 

The  expressions  for  the  stress  and  displacement,  which  have  been 
found  in  this  section, will  allow  the  energy  lost  at  the  moving  crack 
tip  to  be  evaluated.  Thus,  the  analysis  now  proceeds  to  an  energy 
balance. 

B.   ENERGY  BALANCE 

An  energy  balance  will  now  be  applied  to  the  entire  specimen.  First 
the  specimen  energies  are  defined  as  follows:  W  is  the  work  input  to 
the  specimen  by  forces  acting  on  the  boundaries;  I  is  the  reversibly 
stored  energy  in  the  viscoelastic  specimen;  D  is  the  energy  which  is 
dissipated  throughout  the  specimen  due  to  the  viscous  action  of  the 
material;  finally,  S  is  the  energy  which  is  dissipated  at  the  crack 
tip  in  the  creation  of  new  surface  area. 

Neglecting  heat  and  inertia  terms,  the  First  Law  of  Thermodynamics 
may  be  written: 

W  =  I  +  D  ■  +  S 
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The  significant  variables  in  the  problem  are  c  and  t,  thus  the  power 
equation  becomes: 


3W 
ac 


c  + 


3W 

at 


al 

al 

3D 

3D 

aS 

aS 

c 

+ 

+■ 

c 

+ 

+ 

c 

+ 

ac 

t 

at. 

c 

ac 

t. 

at 

c. 

ac 

t 

at 

Rearranging  terms  the  following  equation  is  obtained. 


3W 
ac 


ac 


3D 

3C 


as 

3C 


c  = 


Jt 


3W  +  _3_j_  +  3D  +  3S 
3t    3t    3t    3t 


(ID 


The  left  hand  side  of  Equation  (11)  may  be  thought  of  as  describing 
the  elastic  behavior  of  the  viscoelastic  specimen  (that  behavior  which 
is  independent  of  time)  while  the  right  hand  side  describes  the 


as 


=  0  since  there 


viscoelastic  (time  dependent)  behavior.  However,  -rr- 

a  U 

is  no  mechanism  to  create  energy  S.  Thus,  for  a  fixed  crack  length, 
the  first  law  becomes  W  =  I  +  D  and  the  right  hand  brackets  of 
Equation  (11)  sums  to  zero.  Also,  viscous  dissipation  is  physically 
an  energy  which  must  occur  over  finite  time,  thus  in  the  left  hand  side 
of  Equation  (11 ) : 


3D 
3c 


0 


Therefore  if  time  is  held  constant,  Equation(ll)  becomes 


3W 

3C 


il' 

3C 


3S 
3C 


=   0 


(12) 


Of  the  energies  which  appear  in  Equation  (12),  W  is  the  most 

directly  evaluated.   It  is  simply  the  energy  produced  by  the  applied 

force,  a  ,  moving  through  the  displacement  at  the  boundary.  The 
a 

remaining  energies  of  Equation  (12)  will  now  be  considered  in  detail. 
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1 .  Reversibly  Stored  Energy 

To  further  clarify  the  definition  of  I  at  any  time  t,  consider 

the  following.  Uniaxial  viscoelastic  behavior  can  be  modeled  by  the 

three-element  mechanical  model  consisting  of  a  Maxwell  element  and  a 

spring  loaded  in  parallel  as  shown  in  Figure  3.  The  short-time  spring 

constant  for  the  model  is  E„  =  E  +  E^  where  E  >>  E  .  The  spring 

g       r  r  J 

constants  of  the  model  correspond  to  the  slopes  in  a  stress-strain 
plot  while  model  force  and  displacement  correspond  respectively  to 
viscoelastic  stress  and  strain. 
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Figure  3.  Stress-Strain  Plot  and  Mechanical  Model  for  Visco- 
elastic, Reversibly  Stored  Energy. 

Consider  the  loaded  viscoelastic  body  (model)  at  any  time  t 

(point  A).  Allow  an  instantaneous  relaxation,  thus  losing  no  energy 

in  viscous  dissipation  in  the  model  dashpot,  and  calculate  the  energy 
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released.  Since  this  is  a  conservative  process,  the  energy  released 

will  be  the  reversibly  stored  energy. 

The  relaxation  occurs  with  the  short-time  modulus  to  a  point 

af  zero  stress  (point  B).  The  energy  released  per  unit  volume  is  the 

area  under  the  path  AB  in  the  stress-strain  plot  or: 

1 :     1  a2(t) 
VOL     2   E 

g 

Comparing  this  to  the  strain  energy  of  a  uniaxially  loaded  elastic 
body: 

U  1      a2 


VOL  2     E 

the  following  correspondence  may  be  made  between  elastic  strain  energy 
and  internal  energy  in  specimens  of  identical  geometry  which  are  loaded 
uniaxially. 

U  ■*  I 
ff  E  +  Lg  (13) 

and  a  ■*■     (t) 

2.  Energy  Dissipated  at  the  Crack  Tip 

The  change  in  energy  S  with  respect  to  length,  at  the  crack  tip, 
required  in  equation  (12),  is  now  evaluated.  Holding  t  =  t  ,  allow  a 
small  propagation  of  a  (see  Figure  4).  Since  this  is  for  constant 
time,  the  behavior  of  the  viscoelastic  material  is  elastic  and  the 
evaluation  technique  used  by  Irwin  [Ref.  9]  for  elastic  materials  may 
be  applied.  Thus,  the  energy  lost  in  this  process  is: 


AS  =  I        f    u     an    dx' 
o  J         nmax   max 
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c(t-) 


Figure  4.  Crack  Tip  for  Incremental  Propagation,  a,  at  t  =  t  . 

Substituting  for  the  maximum  values  of  displacement  and  stress,  this 
becomes: 


AS  = 


1 


f   un  <x'  "  a 


tc)  an  (x\  tc)  dx' 


or  in  inertial  coordinates: 

i        /*c(tc)  +  a 

AS  =  7  J   un(x'c(tc}  "  a>  tc)  an  (x  '  c(tc}>  tc)  dx   (14) 

c(tc) 

Equation  (14)  requires  that  the  displacement  of  Equation  (9)  be 
translated  by  x  =  a.  Note  that  the  lower  limit  on  the  integral  of 
Equation(lO),  which  is  a  function  of  position,  thus  changes  to: 


*    +  X  -  C  (t)  -  a   _ 

m 


1  "  t 

m 


(15) 
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With  this  change  in  limits,  the  displacement,  Equation  (9),  and  the 
stress,  Equation  (4),  are  substituted  into  Equation  (14)  to  give: 


AS  = 


w 


c(tc) 


rc(t  )+a     .t. 


t- 


K^t)  /c(i)  -  X  +a 


'm 


dxdx 


Dcr(tc-x) 

3 

3.T 

/  x-c(tc) 

(16) 


Consider  the  compliance  in  Equation  (16).  At  the  upper  limit 
in  time  it  becomes  D  (0)  while  at  the  lower  limit  it  is: 


D 


c(tc)  -  x  +  a 


cr 


4m 


which  is  a  function  of  x.  Now  evaluating  at  the  upper  limit  in  x,  it 
again  is  D  (0)  while  at  the  lower  limit  is  becomes  D  (a/c  ).  Thus 
since  a  is  small,  the  compliance  may  be  approximated  by  the  constant 
value,  Dcr(^)-  The  integrand  in  time  is  now  a  perfect  derivative  and 
Equation  (16)  becomes: 

w 


C(t  )+a 


A3  = 


ocr(o) 


c(tc) 


^? 


m         m 


-  x  +  a 


^{tQ)   /c(tc)  -  x  +  a 


dx 


(17) 


Consider  the  integrand  of  Equation  (17)  which  is  evaluated  at 

any  x  for  which  c(t  )  *  x  *  c(t  )  +  a.  Choosing  x,  c(t  -  fr)   is  seen 

to  be  equal  to  the  crack  length  at  the  time  it  passed  x  minus  the  length 

of  time  it  took  to  propagate  through  a  distance  a  at  its  mean  velocity, 

c  .  In  other  words,  c(t  -  ^— )  is  equal  to  x  -  a  and  the  second  term  in 

m 
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the  brackets  vanishes.     Now  the  time  dependence  in  this  energy  has  been 
evaluated  and  Equation   (17)  may  be  conveniently  written  in  local 


coordinates  as 


AS  =  -!— ^-  Dcr(0) 
it 


ta, 


which  when  integrated  becomes: 

K?(t  ) 
AS  =-i— £-  Dcr(0)  a 


Since  a  is  small,   the  change  in  S  with  respect  to  c  is: 

3S  AS  §#   Dcr'°)  ..-, 

lc    '    7"    =    ■ (18) 

where  the  subscript  on  t  has  been  dropped  since  t  was  a  typical  time. 
3.  Specific  Specimen  Geometry 

At  this  point,  the  energies  required  in  Equation  (12)  have  been 
defined  or  evaluated.  However,  since  W  and  I  depend  on  the  particular 
geometry  and  loading,  the  evaluation  for  these  energies  has  been  only 
general.  A  specific  case  is  shown  in  Figure  5  where  the  applied  load, 
aa(X,t),  is: 

oa(x,  1  b,t)  =  iaa  for  x  >  c(t)  J  (19) 

(0   for  x  <  c(t)) 

The  strains  in  the  specimen  may  be  thought  of  as  a  superposition 
of  the  uniaxial  strain  for  x  ;>  c(t)  and  the  strain  produced  by  the 
stress  field  in  the  vicinity  of  the  crack  tip.  Since  the  uniaxial 
displacement  of  the  boundary  edge  is  e(t)b,  choosing  b  sufficiently 
large  will  make  this  displacement  of  higher  order  than  that  produced  on 
the  boundary  by  the  stress  field  at  the  crack  tip.  Thus  holding 
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Figure  5.  Specimen  Geometry  and  Loading. 

time  constant  and  allowing  an  incremental  propagation  of  dc,  while 
removing  load  from  the  boundary  as  specified  by  Equation  (19),  produces 
negligible  displacement,  hence  negligible  work  at  the  boundary. 
Therefore: 


3W 
3C 


(20) 


For  this  specimen,  if  the  uncracked  length  1  -  c(t)  is  maintained 
large  enough  so  that  a  strip  at  the  right  hand  edge  is  under  uniaxial 
stress,  that  is,  unaffected  by  the  crack  tip  stress  field,  Ref.  11  gives 
the  change  in  the  elastic  strain  energy  as: 


au 

3C 


-rb 


(21) 


Thus  by  applying  Equation  (13) 


3c 


(22) 
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Substituting  Equations   (18),   (20)   and   (22)   into  Equation   (12),  noting 
that  now  surface  energy  is   released  from  two  surfaces  and  that  E     = 
1/Dcr(0),  gives: 


a 


(23) 


Note  that  since  a     is  constant,  K-,  is  a  constant.  This  result 
a  l 

is  identical  to  that  which  is  obtained  in  the  elastic  case. 


C.  ADDITIONAL  RESULTS 

Lindsey  [Ref.  11],  applying  Graham's  "extended  correspondence 

principle,"  transformed  the  elastic  expressions  for  stress  intensity 

factors  directly  into  the  viscoelastic  case.  His  results  agree  with 

Equation  (23)  and  include  the  additional  case  of  a  strip  with 

displacement  boundary  conditions  under  a  constant  displacement,  u^ 

a 

(see  Figure  6).      (Uncracked  specimens  of  this  geometry  are  commonly 
used  in  the  propellant  industry  and  are  referred  to  as   "biaxial   strips" 
since  they  are  in  a  state  of  biaxial    stress.)     Lindsey's  solution   is: 


K 


■  °(t)VT" 


(24) 


eft) 


<r(t) 


IX 


Figure  6.     Cracked  Biaxial    Strip. 
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D.     TIME-TEMPERATURE  SHIFT 

The  expressions   for  K-.(t),  Equations   (23)  and  (24)  are  for  a 
constant  reference  temperature,  T  .     Since  the  viscoelastic  material 
was  assumed  to  be  thermorheologically  simple,     Ref.   12     shows  that  the 
modulus,  hence  the  stress,     may  be  shifted  from  any  temperature,  T, 
to  the  value  at  T    by  multiplying  it  by  the  factor  T    p  /T  p,  where 
p  is  the  material   density.     The  time  scale  is  corrected  by  dividing  it 
by  a  time-temperature  shift  factor,  a-r.     Since  p  -  p   ,  the  expression 
for  the  stress   intensity  factor  become: 

Kl(l7)=    r    °a(y  ^~  stress  B.C.  (25) 

Ki  (if)=  r°(-yVF  dis"-B-c-  <26> 

These  expressions  can  now  be  used  to  guide  experimental  measurements 
of  K-,(t/aT).  For,  in  light  of  the  previous  discussion  concerning  the 
relationship  between  K-,  and  c,  laboratory  measurements  of  these  two 
parameters  in  fracture  tests  will  constitute  a  viscoelastic  fracture 
characterization  of  the  material. 
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III.   EXPERIMENTAL  RESULTS 

Experimental  testing  was  conducted  on  a  CTPB  solid  propellant. 
This  filled  material  consisted  of  solids  (oxidizer)  comprising 
approximately  85%   of  the  total  volume,  bonded  in  a  rubbery  (fuel) 
matrix. 


A.   CONSTANT  STRESS  CASE 

The  specimen  configuration  to  achieve  the  stress  boundary  conditions 
of  Figure  5  is  shown  in  Figure  7.  The  specimen  was  cut  from  samples 
designed  for  use  in  biaxial  tests  with  a  gage  length  of  1.0  inch  and 
thickness  of  0.1  inch.  The  wooden  tabs  were  segmented  with  saw  cuts 
which  stopped  in  the  adhesive  bond  just  short  of  the  propellant. 
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Specimen  Configuration  for  Stress  Boundary  Conditions 
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(To  penetrate  through  the  bond  into  the  propel! ant  would  have  created 

a  stress   riser  and  a  point  of  probable  crack  initiation.)     The  width, 

w,   ideally  would  be  very  small,  allowing  shear-free  edges;  however, 

practical   considerations  forced  a  width  of  w  =0.25  inch. 

Initial   experiments,  with  the  right-hand  edge  of  the  propel  1  ant 

flush  with   the  wooden  tab  edge,  showed  that  soon  after  loading  a  fast 

running  crack  would  initiate  at  a  local   defect  along  the  right  edge  of 

even  a  carefully  cut  specimen.     Adding  a  1/8  inch  unloaded  length  reduced 

the  stress  along  this  edge  and  eliminated  the  problem.     To  compensate 

for  this  departure  from  uniform  stress  at  the  right  hand  edge  and  yet 

still  meet  the  restriction  leading  to  Equation  (21),  c(t)         was  limited 

3  ^  v     /,     \   /max 

to  2.0  inches. 

The  thickness  of  the  saw  cuts  in  the  wooden  tabs  was  large  enough  to 

allow  unrestrained  e    ,     or  d/w  >  ev         .     The  upper  segments  were 

x  x  max 

individually  suspended  and  the  lower  segments  were  individually  loaded 
by  weights  W.  An  initial  crack  of  one  inch  was  cut  with  a  razor  blade. 

A  specimen  was  tested  by  applying  all  loads  at  t  =  0.  The  length 
of  the  propagating  crack,  c(t),  was  measured  with  an  optical  comparator. 
As  the  crack  passed  a  segment,  the  weight  on  that  segment  was  removed 
in  proportion  to  the  width  of  the  segment  traversed  by  the  crack.  The 
tests  were  conducted  in  an  environmental  chamber  with  the  relative 
humidity  less  than  30%,  at  six  temperatures  ranging  from  15°F  to  85°F. 

The  stress  intensity  factor  for  each  test  was  calculated  by  using 
Equation  (25),  assuming  a-r  =  1 ,  and  plotted  on  a  graph  of  stress 
intensity  factor  versus  crack  velocity,  for  each  specific  testing 
temperature  (see  Figures  8  through  13).  In  this  way  the  shift  due  to 
temperature,  in  the  ordinate,  was  accounted  for.  However,  complete 
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shifting  of  the  function  required  that  a-r  be  known.  This  was  obtained 
graphically  by  shifting  the  curves  for  various  temperatures  along  the 
abscissa,  relative  to  the  curve  at  T  ,  so  that  a  continuous  function 
was  obtained.  The  resulting  characterization  curve  is  shown  in  Figure 
14,  where  TQ  =  70°F. 

B.  CONSTANT  DISPLACEMENT  CASE 

Additional  tests  were  conducted  on  cracked  biaxial  strip  specimens. 
Though  the  ends  were  free,  rather  than  restrained  as  shown  in  Figure  6, 
the  long  length  relative  to  the  width  allowed  the  actual  specimen  to 
closely  approximate  the  theoretical  one.  The  strips  were  six  inches 
in  length,  one  inch  in  width  (2b)  and  1/10  inch  thick.  Redwood  tabs 
bonded  to  the  horizontal  edges  allowed  for  gripping.  An  initial  crack 
of  0.75  inch  was  cut,  as  suggested  by  Mueller  [Ref.  7].  As  before, 
the  length  of  the  propagating  crack  was  measured  as  a  function  of  time 
with  an  optical  comparator. 

The  testing  was  performed  with  the  apparatus  which  was  used  for 

similar  tests  on  solithane,  as  described  in  Ref.  13.  Briefly  this 

consisted  of  an  Instron  Universal  Testing  machine  equiped  with  special 

jav/s  to  clamp  the  specimen  and  apply  displacement  u  while  maintaining 

a 

the  upper  and  lower  boundaries  parallel.  The  displacement  was  applied 
at  rates  ranging  from  one  inch/minute  for  low  temperatures  to  five 
inches/minute  at  the  high  temperatures.  A  conditioning  chamber  allowed 
the  tests  to  be  conducted  at  the  same  specific  temperatures  as  in  the 
previous  case. 

Using  Equation  (26),  K,  was  calculated  and  plotted  in  Figures  8 
through  13.  Shifting  was  performed  as  in  the  previous  case,  with  TQ  = 
70°F,  and  the  master  characterization  curve  plotted  in  Figure  14. 
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IV.  DISCUSSION  OF  RESULTS 

In  the  stress  boundary  condition  case,  the  crack  tip  stress 
intensity  is  constant.  This  constant  stress  intensity  should  produce 
a  constant  velocity  crack,  a  fact  which  was  verified  in  these  experi- 
ments. Therefore,  each  sample  tested,  for  a  given  loading  and 
temperature,  produced  a  unique  point  on  a  plot  of  K-,(c). 

However,  for  the  displacement  boundary  condition  case,  since  the 
stress  relaxes,  the  crack  will  decelerate  until  E"  (rubbery  modulus) 
is  reached.  For  the  propellant  tested,  this  relaxation  time  was  on  the 
order  of  a  year.  Therefore,  for  all  experiments  of  this  case,  the 
crack  was  constantly  decelerating.  Thus  K-,(c)  was  a  continuous  function 
and  in  this  study  was  shifted  more  readily  than  the  stress  data,  since 
for  the  later  case,  limited  samples  precluded  enough  testing  to  obtain 
a  well  defined  function  for  each  temperature.  For  this  reason  the 
shift  factors  obtained  in  the  displacement  case  were  used  in  the  stress 
case  to  shift  the  data  to  the  master  curve  shown  in  Figure  14. 

The  theory  of  Section  II  considered  a  purely  viscoelastic  material. 
However,  the  solid  propellant  was  a  filled  material,  which  when  loaded, 
exhibits  the  additional  mechanism  of  dewetting.  (Dewetting  is  the 
breaking  of  bonds  between  solid  filler  particles  and  the  rubbery 
matrix.)  Wood  [Ref.  14]  has  shown  that  propellant  dewetting  is  not  a 
rate  process,  rather  it  occurs  upon  initial  application  of  stress. 
This  behavior  has  a  significant  bearing  on  the  experimental  results 
obtained  in  this  study  as  will  now  be  shown. 

In  both  cases  tested,  the  measured  load  parameter  was  stress,  thus 
automatically  accounting  for  the  modulus  relaxation  in  the  displacement 
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case.  However,  the  two  stress  histories  were  quite  different  as  shown 
in  Figure  15.  In  the  stress  case,  stress  was  constant  while  in  the 

displacement  case  there  was  an  initial   extremely  high  stress  which 

3 
continually  relaxed. 


-*  t 


Constant  Stress 


Constant  Displacement 


Figure  15.  Stress  Histories. 

Thus  a  sample  tested  in  the  stress  case  exhibited  minimum  dewetting 

because  of  its  history  of  constant  stress.  However,  in  the  displacement 

case,  for  a  corresponding  stress  at  t  >  lOt  (where  t  is  the  time 

a       a 

required  to  apply  the  load)  much  more  dewetting  had  occurred  due  to  the 
initial  high  stress.  Consequently  the  crack  propagated  through  a 
material  which  had  been  significantly  altered  from  that  in  the  stress 
case.  The  propagating  crack  had  fewer  bonds  to  break  per  unit  area, 
compared  with  the  stress  case,  and  accordingly  its  velocity  was  as  much 


The  theory  of  this  paper  does  not  account  for  any  transients, 
rather  considers  only  step  loading.  However,  the  application  of  the 
loading  in  the  displacement  case  required  finite  time.  Therefore,  data 
taken  for  this  case  were  taken  at  t  >  10ta  assuming  that  at  this  time, 
any  effects  of  the  loading  method  had  died  out. 
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as  two  orders  of  magnitude  higher.  This  produced  significantly 
different  characterization  curves  for  the  two  cases  as  shown  in 
Figure  14. 

Since  dewetting  produces  microscopic  voids  in  the  material,  stiff- 
ness is  reduced  as  dewetting  increases.  Using  this  fact,  the  increased 
dewetting  of  the  displacement  case  was  verified  in  preliminary  tests  at 
15°F.  In  these  tests,  for  identical  stress  and  time,  the  modulus  in 
the  displacement  case  was  one-half  that  of  the  stress  case. 

In  an  attempt  to  reduce  the  initial  high  stress  peak  in  the  dis- 
placement case,  thus  reducing  the  dewetting,  additional  preliminary 
tests  at  15°F  were  conducted.  The  maximum  stress  was  limited  to  that 
which  produced  stress  intensity  factors  comparable  with  those  of 
the  stress  case  for  the  same  temperature.  Equivalent  moduli  showed 
that  the  amount  of  dewetting  was  in  fact  equivalent.  However,  as  soon 
as  the  displacement  rate  was  stopped,  the  stress  relaxed  due  to  E  -.(t). 
Thus  even  if  the  crack  propagated  at  the  same  velocity  as  in  the  stress 
case  (approximately  0.1  inch  per  hour)  relaxation  of  the  driving 
stress  decelerated  the  crack  to  zero  velocity  over  distances  too 
small  to  be  accurately  measured.  Thus,  no  additional  characterization 
data  were  obtainable  with  this  loading  method. 

The  accuracy  of  the  measured  velocity  in  the  stress  case  was  con- 
sidered superior  to  that  of  the  displacement  case,  especially  for  low 
velocities.  For  in  the  displacement  case,  c(t)  had  to  be  measured  over 
relatively  short  periods  of  time  to  obtain  an  accurate  crack  propagation 
function.   (Velocity  was  obtained  by  graphically  measuring  the  slope  of 
this  function.)  As  C(t)  became  small,  Ac(t)  for  the  time  interval 
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mentioned  above  also  became  small.  Thus,  because  the  experimental  crack 
propagation  was  irregular  due  to  inhomogenei ties  in  the  propellant,  the 
actual  velocity  was  obscured  as  those  irregularities  approached  the 
magnitude  of  Ac(t).  However,  in  the  stress  case,  c(t)  was  a  constant 
and  Ac(t)  could  be  measured  over  large  time  intervals. 

The  scatter  in  the  characterization  curves  is  due  to  the  inaccur- 
acies mentioned  above  for  the  displacement  case  and  also  to  the 
inhomogeneous  nature  of  a  highly  filled  viscoelastic  material  such  as 
propellant.  By  comparison,  similar  characterization  of  solithane,  an 
unfilled  and  relatively  more  homogeneous  material,  by  Francis  et.  al . 
[Ref.  13]  showed  significantly  less  scatter. 

Referring  to  the  characterization  curves,  in  both  cases  there  appears 
to  be  a  K,  critical  below  which  a  crack  will  not  propagate.  Again  for 
both  cases,  especially  apparent  in  the  stress  case,  there  appears  to  be 
a  maximum  velocity  of  propagation;  however  these  maximum  velocities 
are  somewhat  misleading.  Preliminary  tests  showed  that  at  very   high 
loads  the  failure  was  general  and  catastropic,  initiated  at  voids  and 
other  stress  risers  throughout  the  specimen  as  well  as  at  the  crack  tip. 


Solithane  at  room  temperature  has  a  relaxation  time  of  less  than 
one  minute.  Thus,  though  tested  for  the  displacement  case,  essentially 
constant  velocity  crack  propagation  in  these  tests  allowed  the  additional 
advantage  of  increased  accuracy  of  the  measured  crack  propagation 
function  for  reasons  mentioned  above. 
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V.  CONCLUSIONS 

The  theory  presented  in  this  paper  permits  fracture  characterization 
of  a  viscoelastic  material  to  be  made.  However,  for  a  filled  visco- 
elastic  material  such  as  a  solid  propellant,  dewetting  behavior  affects 
the  crack  velocity.  (The  stress  boundary  condition  case  is  the  limit 
case  of  minimum  dewetting.)  Additional  theory  is  necessary  to  account 
for  dewetting.  This  theory  would  then  be  expected  to  provide  for  a 
superposition  of  the  displacement  characterization  curve  upon  the  stress 
characterization  curve  (See  Figure  14). 
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